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ABSTRACT 
The  primary  problem  in  the  calculation  of  the  differential 
scattering  cross  section  for  hydrogen  molecular  ions  is  the  derivation 
of  an  approximate  expression  for  the  probability  that  the  hydrogen 
molecular  ion  will  be  oriented,  prior  to  impact,  in  such  a  way  as  to 
avoid  disruption  during  scattering.   Using  a  simple  model,  an  approxi- 
mate expression  for  the  orientation  probability  is  obtained. 
Differential  cross  sections  for  the  scattering  of  hydrogen  molecular 
ions  from  carbon  and  krypton  are  calculated  using  classical  scattering 
theory  and  the  derived  orientation  probability.   Impact  energies  of 
100  to  1000  electron  volts  are  considered.   In  addition  differential 
cross  sections  for  protons  scattered  from  carbon  and  krypton  are 
presented. 
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1 .   INTRODUCTION 

Bates,  Cook,  and  Smith,  have  presented  a  classical  theory  of 
ion-molecular  rearrangement  collisions  at  high  impact  energies    . 
Using  proton-methane  collisions  as  an  example  a  brief  outline  of  the 
theory  follows. 

In  Figure  1  m.  is  the  impacting  proton,  mg  is  one  of  the 

hydrogen  atoms  of  the  methane  molecule  and  M,  is  the  carbon  atom  of 

the  methane  molecule.   In  this  collision  m.  =  mg  «  M,  .  For  capture 

/     e 
to  occur  V.  fa   V2  .  From  these  two  conditions  and  the  conservation  of 

of  linear  momentum  and  energy  it  can  be  shown  that  Q2   **  90°  and 
O^  «s  45  .  Thus,  for  this  case,  capture  of  the  hydrogen  atom  of  the 
methane  by  the  incident  proton  to  form  hydrogen  molecular  ions  may 
occur  at  an  angle  of  approximately  45  .   Bates,  Cook  and  Smith  have 
calculated  a  cross  section  for  the  production  of  hydrogen  molecular 
ions  by  this  process 

An  experimental  difficulty  may  arise  in  a  measurement  of  the 
cross  section  for  the  production  of  Hg  by  this  double  scattering 
process  if  the  incident  proton  beam  contains  some  significant  number 
of  Hg .   If  the  cross  section  for  scattering  of  Hg  from  the  carbon 
atom  is  much  larger  thm  the  cross  section  for  H2  production  by  capture, 
then  a  measurement  of  the  production  cross  section  would  be  impos- 
sible unless  a  pure  source  of  protons  were  used.   In  order  to 
determine  the  relative  magnitudes  of  these  two  differential  scatter- 
ing cross  sections  a  calculation  of  the  differential  cross  section 
of  Hg  ions  scattered  from  heavy  atoms  is  presented. 


1>J 


■■> 


•-.» 


■ 


CD 

CD 
CJ 

O 
U 

a 

cu 
u 

■P 

Pu 
CO 
O 

a> 


o 

CJ 

o 

4J 

nJ 
n 
■u 

CD 

3 


>S 


■u 


o 


cu 

M 

00 


10 


2.   THE  MODEL 

The  hydrogen  molecular  ion  is  linear  with  an  inter  atomic  spacing 
of  approximately  two  angstroms.   For  randomly  oriented  impacting  Hg 
ions  this  implies  a  possible  difference  in  impact  parameter  for  the 
individual  atoms  of  as  much  as  two  angstroms.   If  the  difference  in 
impact  parameters  is  great  enough,  the  velocities  of  the  two  atoms 
after  scattering  may  differ  by  enough  to  cause  their  relative  energy 
to  exceed  the  binding  energy  of  the  molecule.   The  molecule  will 
then  be  disrupted  and  can  not  be  considered  an  H2  scattering  event. 

Even  at  zero  difference  in  impact  parameter  disruption  of  the 
molecule  may  occur  because  of  the  following  consideration.   In  a 
spherically  symmetric  potential,  such  as  that  used  in  this  calcu- 
lation, each  particle  remains  in  a  plane,  defined  by  the  scattering 
axis  and  a  line  from  the  scattering  axis  to  the  particle,  during 
and  following  the  collision.  If  the  angle  between  the  planes  is 
great  enough  the  difference  in  the  velocities  of  the  two  atoms, 
after  scattering,  may  be  such  that  the  molecule  is  disrupted. 

Figure  2  shows  the  relationship  between  the  scattering  planes 
and  the  scattering  axis.  £<i>  is  the  angle  between  the  scattering 
planes. 

The  condition  that  the  molecule  remains  intact  may  be  stated  as 
D  >  h  M(CV)2  where  D  is  the  binding  energy,  £V  is  the  relative 
velocity  of  the  two  particles  and  %  £t(£V)2  gives  the  relative  energy 
of  the  two  particles.   Using  the  symbol  AV  for  the  maximum  allowable 
£V  the  limiting  condition  for  nondisruption  becomes: 
D  =  k  y,  (AV)2  . 
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Thus  both  the  difference  in  impact  parameter  £p  and  the  angle  between 
the  scattering  planes  must  be  such  that  the  above  condition  is 
satisfied.   This  limitation  on  the  orientation  of  the  impacting 
molecule  leads  to  an  orientation  probability  that  must  be  considered 
in  the  calculation  of  the  differential  scattering  cross  section  for 
hydrogen  molecular  ions. 

For  high  impact  energies  and  small  impact  parameters  each  atom 
of  the  He  may  be  treated  as  a  proton.   That  is,  the  effective  charge 
may  be  taken  as  one  for  each  particle.  The  calculation  of  the  differ- 
ential cross  section  for  H^  ions  may  then  be  approached  in  the  follow- 
ing way.   For  a  given  energy  and  impact  parameter  calculate  the 
scattering  angle,  and  the  differential  cross  section  for  that  angle, 
for  one  proton  whose  velocity  is  V0 ,  the  velocity  of  the  Hg .  Then 
determine  the  probability  that  the  "second"  particle  is  oriented  with 
respect  to  the  "first"  particle  in  such  a  way  that  its  scattering 
angle  is  within  A  0  of  Q   and  its  scattering  plane  makes  an  angle  of 
less  than  A  $  with  the  scattering  plane  of  the  first  particle.  The 
differential  cross-section  for  H2  scattering  will  then  be  given  by: 
Ou*(B)   =  CJut(8)  P,  where  P  is  the  orientation  probability.  This  is 
then  done  over  a  range  of  impact  parameters  and  for  several  energies. 
3.   CALCULATION  of  <7Hf(8) 

The  calculation  of  (7^(8)  for  the  "first"  particle  was  done  using 

the  method  suggested  by  Bates,  Cook  and  Smith   .   The  potential  used 

for  the  calculation  was  a  shielded  potential  of  the  form, 

V(R)  =  —  (Me*"+  Ne'  +  Le'  )  due  to  Byatt   .   Calculations  were  made 
R 

for  both  carbon  and  krypton.   Byatt  lists  several  potentials  in 
reference  2,  all  of  which  may  be  conveniently  used  in  the  computer 
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program  in  appendix  2.   The  cross  section  was  calculated  fox  impact 
parameters  from  0  -  5  a0  in  increments  of  0.05ao  (a0  is  the  radius  of 
the  Bohr  orbit)  and  for  energies  from  lOOev  to  lOOOev  in  steps  of 
lOOev.   The  details  of  the  calculation  are  shown  inappendix  1. 
4.   CALCULATION  OF  THE  ORIENTATION  PROBABILITY. 

The  condition  AV  =/4D  \i   may,  for  convenience,  be  put  in  terms 
of  the  energy  of  the  %.      Bj  =  %(2MH)VJ  =  M^  or  M„  =  Ja  then 

0 

AV  =  /4D  \  2  V0  .   This  limiting  condition  implies  a  relationship 

hi 

between  £Q  and  £$>  such  that  AV=  F(G6,  G$).  £9  is  related  to  £p  by 
QQ   =  pCp/CH*(8)  SINQ  in  the  center  of  mass  coordinate  system. 
Therefore  AV  =  G(£<E>,Gp)  gives  a  relationship  between  Q$  and  Qp.      Since 
AV  is  a  constant  for  a  given  energy  and  impact  parameter,  there  exists 
some  relationship  between  £p  and  C&  such  that  £<i>  =  f  (Cp).   C^  =  f  (Cp) 
defines  a  curve  in  a  plane  perpendicular  to  the  scattering  axis.   The 
second  particle  must  pass  through  the  area  on  the  plane  bounded  by 
this  curve  if  disruption  is  to  be  avoided. 

The  randomly  oriented  H^  can  be  considered  as  defining  a  sphere 
the  center  of  which  is  the  "first"  atom  and  whose  radius  is  X,  the 
inter  atomic  spacing. 

In  Figure  3,  A  is  the  area  on  the  surface  of  the  hemisphere  con- 
tained within  the  surface  defined  by  Cp  =  f  (£$)"•   T^6  ratio  of  this 
area  to  the  area  of  a  hemisphere  of  radius  I   gives  the  fraction  of 
the  Hg  ions  which  are  oriented  in  such  a  way  as  to  avoid  disruptions. 
Thus  the  orientation  probability,  P  is  given  by:   P  =  A/lfTl2  .   No  exact 
solution  has  been  obtained  for  the  orientation  probability. 

As  an  approximation  to  the  orientation  probability  the  following 
simplifications  were  made. 
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Assuming  the  target  atom  to  be  infinitely  massive,  there  will  be 
no  change  in  the  magnitude  of  the  velocity  of  either  atom  of  the 
hydrogen  molecular  ion.   Since  they  then  have  the  same  initial  speed 
V0  ,  \i  =  V2  =  V0  .   When  Q$  =   0  the  particles  are  scattered  in  the  same 
plane  as  shown  in  Figure  4a.   Then,  assuming  A9  is  small 

l^-V  .   When  C9  =  0  (w^icfe  implies  £p  =  0) 

the  relationship  between  the  individual  atoms  is  as  shown  in  Figure  4b. 

^-\w  or  AJ>=/*2_W  1  \  these  approximations 

then  give  independent  maxima  for  £<£>  and  QQ .   £8  is  related  to  /p  by  z. 

Ap  -  V(e)  SINe  &     or  Ao  =  <V(0)  SIN0i 


m 


The  above  restriction  on  Ap  implies  that  the  second  patticle  must 
be  contained  between  two  coaxial  cylinders  whose  axis  is  the  scattering 
axis  and  whose  radii  are  p  -  Ap  and  p  +  Ap  respectively.   The  condition 


s)2pW-) 


puts  a  limit  on  the  angular  separation  of  the  scattering  planes.   The 
limits  on  A<£  and  Ap  then  confine  the  second  particle  to  lie  on  the 
surface  of  a  hemisphere  of  radius  i   contained  between  the  two  scattering 
planes  and  the  two  cylindrical  surfaces.   Figure  5  illustrates. 

This  area,  on  the  surface  of  the  hemisphere,  divided  by  the  area 
of  the  hemisphere  then  gives  an  approximation  to  the  orientation 
probability. 

As  a  further  approximation,  this  area  was  assumed  small  and  taken 
to  be  the  plane  area  given  by  A«  p(2A<£)  (2Ap);  which  then  gives,  as 
an  estimate  of  the  orientation  probability*  P  «  2pA^p    t 
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5.   RESULTS  AND  BISCUSSION 

Everhart,  et  al,  have  calculated  the  differential  cross  section 

[3] 
for  scattering  from  a  screened  Coulomb  potential   .   The  results  are 

presented  in  terms  of  a,  the  screening  length,  given  by: 

a  =$»+z?  1  h 

and  b,  the  collision  diameter,  given  by: 

b  =  2l  Zg  e 

where  /I  =  the  reduced  mass  of  the  system. 

V  =  the  relative  velocity  of  projectile  and  target. 

As  a  check  on  the  calculation  of  the  cross  section  for  proton 
scattering  the  computer  program  used  for  the  calculation  was  run  using 
the  screened  Coulomb  potential  of  carbon  and  a  proton  impact  energy 
such  that  the  ration  b/a  was  equal  to  1.   The  results  of  this  calcula- 
tion were  within  2.34%  of  those  of  Everhart  et  al. 

The  actual  calculations  were  made  using  the  Byatt  potentials  of 
carbon  and  krypton  with  proton  impact  energys  from  50  to  500  electron 
volts.   The  results  are  presented  in  the  graphs  of  Figures  6  and  7. 

The  cross  section  for  the  scattering  of  the  hydrogen  molecular 
ion  for  impact  energies  of  100  to  1000  electrons  volts  (corresponding 
to  proton  energies  of  50  to  500  electron  volts)  then  is  the  cross 
section  for  proton  scattering  multiplied  by  the  orientation  probability, 

These  cross  sections  are  presented  in  the  graphs  of  Figures  8  and  0. 

Under  the  assumption  that  the  speeds  of  the  two  particles  are  the 
same  following  scattering,  the  £Q   used  in  the  orientation  probability 
approximation  is  due  only  to  the  difference  in  the  directions  of  the 
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Figure  6.   Graph  of  the  logarithm  of  the  differential  scattering  cross 
section  versus  scattering  angle  for  protons  on  carbon. 
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Figure  7.   Graph  of  the  logarithm  of  the  differential  scattering  cross 
section  versus  scattering  angle  for  protons  on  krypton. 
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Figure  8.   Graph  of  the  logarithm  of  the  differential  scattering  cross 
section  versus  scattering  angle  for  hydrogen  molecular  ions  on  carbon. 
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Figure  9.   Graph  of  the  logarithm  of  the  differential  scattering  cross 
section  versus  scattering  angle  for  hydrogen  molecular  ions  on  krypton, 
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velocity  vectors.   The  speed  of  a  particle  scattered  through  an  angle 
Q  is  given  by: 

1-  2MtMH     (1-  cosQ) 

L   (vf 


V-  v0 


where  V  =  the  speed  of  a  projectile  particle  scattered  through 

an  angle  0  in  the  CM.  Coordinate  system, 

M  =  Mass  of  the  projectile  particle, 

M  =  Mass  of  the  target  particle 
and  V0sr  initial  speed  of  the  projectile  particle. 

If  the  speeds  of  particles  scattered  through  angles  9  and  Q  +  QQ 
differ  by  enough  that  CVS,  resulting  from  this  difference,  is  greater 
than  AV,  the  maximum  for  non  disruption,  then  the  Hg  will  be  disrupted 
due  to  the  speed  difference  alone.   If  AQ5»A9,  then  AQ  (due  to  the 
difference  in  direction  only)  will  be  the  more  stringent  condition  and 
the  error  introduced  by  ignoring  the  difference  in  speeds  may  be  assumed 

small. 

2M  M 
The  factor   t  n   ,  for  scattering  from  carbon,  is  approximately 

.   The  factor  (1-  cos6)  is  always  less  than  one  for  scattering 
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angles  between  zero  and  90  degrees.   Thus  the  expression  for  V  above 

can  be  expanded  in  a  series.   Taking  only  the  first  two  terms  gives, 


V-  V0 


L  MA         (i-  cose) 

(Mt  +  M/ 


] 


M.    SI. 
then  CV     =  -V0  pVjiTg-     sin0  £0 


when  CV5    =  AV  and  C6S     -  A4     then    Ia6c  I    =1^  +^)        /^V_\ 

$  VtsinetM 
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but  AV  = 


(^—  ert- 


Therefore  AV  =  /4D  \  2   ,  which  is  just  the  AQ  due  to  the  differ- 
ence  in  direction  of  the  velocity  vectors. 

Therefore  A9  =  V_5 — Li_  .   For  carbon,  then  A95  ^  14  A9. 
HJl   sin9 

This  implies  that  the  Hg  will  be  disrupted  due  to  the  difference  in  the 
directions  of  the  velocities  of  the  two  particles  at  a  A9  which  is  much 
less  than  the  AQ.  required  to  cause  a  disruption  due  to  a  speed  differ- 
ence only. 

Because  there  is  some  difference  in  the  speeds  of  the  two  particles 
the  relative  energy  of  the  particles  is  somewhat  greater  than  that 
predicted  with  the  speeds  assumed  equal.   Thus  disruption  will  occur 
more  often  than  predicted  by  the  model  used  for  the  calculation. 
Therefore,  the  assumption  of  equal  speeds  tends  to  make  the  probabil- 
ity factor  somewhat  too  large. 

The  use  of  separate  maxima  for  £*p   and  A$  introduces  an  error  in 
that  the  final  probability  expression  contains  the  product  of  the 
maximum  values,  each  of  which  was  derived  under  the  assumption  that 
the  other  was  zero.   This  error  tends  to  make  the  probability  factpr 
too  large. 

In  view  of  the  small  values  of  P,  (representative  values  are 
listed  in  Table  1.)  the  assumption  that  the  area  on  the  surface  of  the 
hemisphere  can  be  approximated  by  a  plane  area  would  seem  to  be 
reasonable  for,  at  least,  the  higher  angles  and  energies. 

The  conclusion  is  that  the  probability  factor  is  in  general  some- 
what too  large  thus  causing  the  calculated  cross  sections  for 
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Energy 
(in  ev) 

6  (Lab) 

Carbon 
P 

Kryp 
6  (Lab .) 

ton 

P 

100 

83.65 

1.95xl0"3 

87.88 

5.08x10 

100 

45.62 

7.72xl0"3 

45.23 

2.22x10 

100 

10.68 

2.10X10'1 

10.32 

6.38x10 

500 

88.88 

9.18xl0"5 

83.45 

2.68x10 

500 

42.35 

7.08xl0"4 

45.39 

1.66x10 

500 

10.37 

2.30xl0"2 

10.10 

7.00x10 

1000 

75.20 

1.72xl0'5 

86.47 

6.04x10 

1000 

47.62 

7.12xl0"5 

46.80 

3.12x10 

1000 

10.35 

7.46xl0"3 

10.13 

2.32x10 

-2 


-4 


-3 


-2 


-5 


-4 


Table  1.   Representative  values  of  the  orientation  probability. 
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hydrogen  molecular  ion  scattering  to  be  too  large. 

In  Figure  10  the  curve  marked  (JH*(45  )  is  the  differential  scatter- 
ing cross  section  for  He  from  carbon  at  a  scattering  angle  of  45  as 
a  result  of  this  calculation.   The  impact  energies  plotted  are  from 

lOOev  to  lOOOev.   The  curve  labeled  aOS+CH^ *  CH3+Hg)  is  the  cross 

section  for  the  capture  of  a  hydrogen  atom  from  the  methane  molecule 
by  an  impacting  proton.   This  curve  is  a  result  of  calculations  based 
on  the  reaction  cross  section  as  derived  by  Bates,  Cook  and  Smith   . 
Figure  10  shows  the  cross  section  for  scattering  of  Hg  from  carbon  at 
45?  to  be  somewhat  more  than  an  order  of  magnitude  greater  than  the 
cross  section  for  the  reaction  H  +CH4 — >CH3+Hg.   As  a  result,  in 
order  to  measure  the  reaction  cross  section  it  would  appear  that  a 
proton  beam  free  of  Hg  will  be  necessary. 
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o 


-f 


LOQ0  Energy  (in  electron  volts) 

Figure  10.   Comparison  of  reaction  and  scattering  cross 
sections . 
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APPENDIX  1 
CALCULATION  OF  &H+ (6) 

The  classical  differential  scattering  cross  section  in  the  center 

of  mass  coordinate  system  is  given  by: 

-1 


where  u=  (-_ dl 


CJCe)  =  e  |>we|g|]""  (1) 

where?  is  the  impact  parameter  and  6  is  the  scattering  angle  in  the 
center  of  mass  coordinate  system,  ©is  related  to(?by:  ©=  7T  —  2€i)      (2) 

rzV{       y"(r)   IjeTP5  m 

in  which  fe  is  the  largest  positive  root  of 

and  s&t   =  the  reduced  mass  of  the  system 
VOO  =  the  scattering  potential. 
Vo  =  the  velocity  of  the  impacting  particle  in  the  lab 

coordinate  system  when  the  target  particles  are  considered 

to  be  at  rest. 

[2] 
The  potential  as  given  by  Byatt,    is 


VCD  =  x'2>3-a  (oe*r  +  e  e'V  f  eer) 


r  (5) 

after  converting  V(f)  to  electron  volts  and  letting  f=RG[0  there  results: 

V(A)  =  i(Aew^  6epR  +cey")  (6) 

Where                    For  Carbon  For  Krypton 

A=                  204.136  328.635 

B=                  -71.856  588.600 

C=                   31.029  63.765 

a=                   -1.700  -1.086 

j3=                  -2.895  =4.961 

y=                                             -8.808  -26.110 

The  quantity  J^V,  can  be  put  in  terms  of  the  energy  of  the  H* 

in  electron  volts  in  the  lab.   coordinate  system. 
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ytC   S  .,  H  T-  in  which  the  mass,  MTlis  the  mass  of  the  target  atom. 
DeflM  ET=   kCTgSlf-)^-  (7) 

(8) 

-  *  "  *%  =  o 


Define     V  =      VMft 
ET 

Then  /?0  is   the  root  of 

T    "    "ft*    =    0  (9) 

The  root  of   this   equation  was    found  using   the  Newton- Raphson  method 


which  is         fi     ■  *„.'-  £CM 


where 


and 


fft)  =    I  -  f  -  * 


fW  =   v^v   t   itf 


(10) 


(11) 


(12) 


in  which  DV  is  defined  by 


then 


6f\ 

The  integral,  U,  was  evaluated  using  the  Gauss-Mahler  Quadrature 


(13) 
(14) 


formula  as  employed  by  Bates,  Cook  and  Smith 


Ll] 


The  result  is 


where 


in  which  the 


To  calculate  JsLzJ  differentiate  equation  2  with  respect  to  $ 
id  Pi 


(15) 
(16) 

(17) 


it 

de 


(18) 
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where       du  .  iy  +  MM* 
df      3?       3/?t,de 


and  the  X.  are  as  defined  previously. 
TT"  is  obtained  by  differentiating 

the  result  is 


in  which 


where  the  X.  are  as  defined  previously, 


(19) 


in  which  »  ..    ,  -  _-| 


(22) 


§r  =  #2  -  v/f-x<  -  •*  ]'[%-*-  **?*]       (24) 
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EXPLANATION  OF  PROGRAM 

STATEMENT  INDEX 

1-5  Places  and  prints  column  headings  for  print  out. 

6-8  Dimension  statements.   Self  explanatory. 

9-23  Do  loop  for  graph  title. 

24-29  Constants  for  the  potential  used.   See  equation 

No.  6,  Appendix  1.   May  take  on  any  values  in 
accordance  with  potential  used. 

30  Sets  the  initial  value  for  the  number  of 
iterations.   See  statement  index  74  and  75. 

31  Sets  the  initial  value  of  the  impact  parameter 
such  that  RH6(1)  =  9.   See  statement  index  60 

32  Sets  the  initial  value  of  theta,  the  scattering 
angle  in  C.  M.  Coordinates. 

33  Sets  the  initial  value  of  the  energy  of  the 
hydrogen  molecular  ion  in  the  LAB.   Coordinate 
system,  EV,  such  that  EV(1)  =  lOOOev.   See 
statement  index  37. 

34,  39,  40,  41,     Give  the  draw  program  the  proper  parameter  to 
42,  43,  44,  45      draw  10  curves  on  one  graph. 

35  Begins  outer  do  loop  for  10  different  energies 
defined  by  statement  index  33  and  37. 

36  Subscript  relationship  to  achieve  desired  energies, 

37  Energy  defining  relationship.   May  be  varied  for 
desired  energies. 

38  Sets  the  initial  value  of  R  for  the  iteration 

to  find  the  root  of  equation  4,  Appendix  2.   May 
be  based  on  b,  the  collision  diameter.   (Pag el  9) 

47-57  Label  program  for  the  labeling  of  curves  with 

appropriate  energy. 

58  Begins  inner  do  loop  for  101  values  of  the  impact 
parameter,  RH0,  as  defined  by  statement  index  31 
and  60.   RH0  is  in  units  of  a^ 

59  Subscript  relationship  to  achieve  desired  values 
of  RH0. 


38 


60  Defining  relationship  for  RHO.  May  be  varied 

as  desired. 

61-63  Evaluates  the  exponential  terms  of  the  potential 

used. 

64  Mass  of  the  projectile  particle  in  a.  m.  u. 

65  Mass  of  the  target  particle  in  a.  m.  u. 

66  Evaluates  the  relative  energy  factor  of  equation 
4,  Appendix  1  in  terms  of  EV. 

67  Evaluates  V  as  defined  by  equation  8,  Appendix  1. 

68  Evaluates  equation  11,  Appendix  1,  for  the 
iteration. 

69  Evaluates  DV  as  defined  by  equation  13,  Appendix 
1,  for  the  iteration. 

70  Evaluates  equation  12,  Appendix  1. 

71  Computes  a  new  value  of  R  for  the  iteration. 
72-73          Compares  the  last  two  values  of  R. 

74  Counts  the  number  of  iterations. 

75  Stops  the  program  if  finding  R  requires  over 
10000  iterations. 

76  Replaces  R  with  the  last  value  computed  by 
statement  index  71. 

77  Repeats  the  iteration  if  statement  index  73  has 
not  been  satisfied. 

78  Continues  the  program  if  statement  index  73  has 
been  satisfied. 

79-81  Insures  that  the  proper  value  of  R  has  been  used 

to  compute  V  and  DV. 

82  Stops  the  program  if  R  is  less  than  RHO  and  exits 

the  inner  do  loop  if  R  =  RHO. 

83-85  Evaluate  the  cosine  functions  of  equation  15, 

Appendix  1. 

86-100  Evaluate  the  g  (X.)  as  defined  by  equation  16, 

Appendix  1. 
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101  Evaluates  the  integral  U  as  defined  by  equation 
15,  Appendix  1. 

102  Computes  theta  the  scattering  angle  in  the  C.  M. 
coordinate  system  for  the  given  energy  and  impact 
parameter  by  equation  2,  Appendix  1. 

103-105  Computes  the  9^'V^P  as  defined  by  equation  21, 

Appendix  1. 

106  Evaluates  3lV3£  as  defined  by  equation  20, 
Appendix  1. 

107  Computes  dRo/j<?  as  defined  by  equation  22, 
Appendix  1. 

108-113  Computes  the  ^»  Vw0  as  defined  by  equations  24, 

Appendix  1. 

114  Computes  S^/^R0  as  defined  by  equation  23 
Appendix  1. 

115  Computes  dU/dp  as  defined  by  equation  19,  Appendix 
1. 

116  Computes  jo^/jo)  as  defined  by  equation  18r  Appendix 

J-  • 

117  Computes  the  differential  cross  section  for 
scattering  of  protons  at  the  previously  computed 
angle  by  equation  1,  Appendix  1. 

118  Computes  A  pas  given  on  page  16. 

119  Prevents  dividing  by  zero  in  statement  index  122 
when  theta  is  180. 

120  A  dummy  statement  to  give  a  value  of  A  $  for 
theta  =  180. 
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121  Skips  statement  index  122  when  theta  =  180. 

122  Evaluates  £$  as  given  on  page  16. 

123  Evaluates  £fi.  as  given  on  page  16. 

124-127  Compute  the  percent  difference  in  velocities 

after  scattering  as  described  on  page  24. 

128-132  Convert  the  scattering  angle,  theta,  to  LAB 

coordinates  and  degrees. 
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133  Gives  the  inter  atomic  spacing  I, 

134  Computes  the  value  of  the  orientation  probability 
as  defined  on  page 

135-136  Limits  the  value  of  P  to  1. 

137  Computes  the  differential  cross  section  for 
Hg  scattering  as  given  on  page  13. 

138  Gives  the  proton  energy  for  printing  when  that 
is  of  interest. 

139  Prevents  the  log.  function  from  having  a  negative 
or  zero  argument  in  statement  index  140. 

140  Computes  the  logrithum  of  the  cross  section  for 
the  draw  program. 

141-142         Limits  the  values  of  the  scattering  angle  printed 
to  those  between  2  and  90  . 

143  Prints  results  for  each  value  of  RHO  assigned. 

144-150         Self  explanatory. 
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